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We theoretically study the quasiparticle interference (QPI) of the surface states in crystalline 
topological insulators which possess mirror symmetry and time-reversal symmetry, by analyzing the 
Fourier transformed local density of states (FT-LDOS), p{c\^uj) around a single static impurity on 
the surface. We show that the symmetries in the QPI pattern are determined by the transformation 
properties of the impurity under mirror reflections and time-reversal. We study the singularities 
in p(q, uj) and show that while the presence of singularities in p(q, uj) depends on the geometric 
features of the iso-energy contour at cj, the absence of certain singularities denotes the scattering 
forbidden by underlying symmetries. We apply the general analysis to the QPI on the (OOl)-surface 
of Pbi-xSua^Te and predict all vanishing singularities in p{c\^uj). The model-independent analysis 
is supported by numerical calculations of p(q,C(;) around a single point charge impurity using an 
effective four-band model. We demonstrate that QPI can also be used to probe the Lifshitz transition 
of the surface states observed in recent ARPES experiment. 



Recent progress in topological materials indicates that 
topological states protected by time-reversal symmetry 
(TRS) may be the first identified class amongst a much 
wider family of topological states protected by discrete 
symmetries'^. Hsieh et at predicted^, via first princi- 
ples calculation, that SnTe is a nontrivial insulator whose 
topological state is protected by mirror symmetries about 
(llO)-planes. The physical manifestation of this nontriv- 
ial topological state is the existence of surface states on 
the (OOl)-plane which contains four Dirac points (two 
along X-T-X and two along Y-T-Y). This prediction has 
been experimentally confirmed^ — through the use of an- 
gle resolved photoemission spectroscopy (ARPES) which 
has successfully identified the predicted surface states in 
Pbi-a^Sua^Te and SnTe. It was also explicitly shownSiii^ 
that as the binding energy increases away from the Dirac 
point binding energy, the surface states undergo a Lif- 
shitz transition at a critical energy Ec where the Fermi 
surface (FS) changes from two separate ellipsoids to two 
concentric ellipsoids. 

While ARPES has established the existence of the non- 
trivial surface states and their dispersion, the informa- 
tion of the wave functions is missing; quasiparticle in- 
terference (QPI) study is a useful way to elucidate this 
information via scanning tunneling microscop}4ii — . In 
this paper, we theoretically examine the symmetry pro- 
tected features in the Fourier transformed local density 
of states (FT-LDOS) around a point impurity, p(q, cj). 
We analyze the general symmetries of p(q, cj), by exam- 
ining the relations among the values of the FT-DOS at 
(Qx.qy), {qx,-qy) and {-Qx.qy), for all types of static 
point impurities. We predict the absence of symmetry 
forbidden singularities in p(q, cj). Generically the sin- 
gularities are given by the stationary points on the iso- 
energy contour^. The absence of any of these singulari- 
ties is determined by the orthogonality of wavefunctions 
at the corresponding stationary points, and is a symme- 



try protected property inherent to the system. We apply 
this result to Pbi-a^Sn^^Te. We design a four-band model 
having all the symmetries of Pbi-a^Sn^^Te and captur- 
ing the overall features of the observed band structure. 
We calculate the FT-LDOS around a point charge im- 
purity. Within the calculated p(q, a;), we can identify 
each singularity with a pair of stationary points on the 
iso-energy contour while symmetry forbidden singulari- 
ties vanish. The QPI pattern changes drastically at the 
Lifshitz transition. 

We begin by considering a bulk 3D insulator whose 
surface is terminated on the (OOl)-plane leaving only 
the in-plane momenta kx and ky as well-defined. Be- 
fore turning to the study of any specific material, we 
want to understand the general features of p(q, u) around 
a point impurity in the resultant 2D system having 
the three symmetries: mirror reflections about xz- and 
yz-pldines and TRS. The three symmetries transform 
the single-particle Hamiltonian, i^(k), are defined as, 
MxzH{kx,ky)M-^ = H{kx,-ky), My,H{kx,ky)M-^ = 
H{-kx,ky), UTH(kx,ky)U:^^ = H^{-kx,-ky), where 

(k) is written in a predefined but arbitrary orbital basis 
and My and Ut {T = KUt) are the unitary matrices 
representing the symmetry operations in the same basis. 
Ut must be antisymmetric for a spinful system. We use 
the exact formulation of T-matrix to calculate p(q, cj) 
associated with any type of single impurities. We may 
express p(q, cj) in terms of the single particle retarded 
Green's function and its complex conjugatei^s^: 



p(q,cj) = 

- 7^ E ^^[^(k' k + q, c^) - G* (k + q, k, uj)] , (1) 
k 

where N is the total number of sites on the surface. The 
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retarded Green's function G(k, k + q, cj) given by 
G(k,k + q,cj) = 

Go(k, uj)S^o + Go(k, cj)T(cj)Go(k + q, cj), (2) 

where Go(k, a;) = (cj — i^(k) + ir])~^^ is the retarded 
Green's function for a homogenous system. The matrix 
T{uj) contains the effect of impurity and is defined as 

T{u;) = V{I-g{uj)V)-\ (3) 

where V is the impurity potential matrix written in the 
orbital basis and g{uj) = Go(k, cj)/7V. A neces- 
sary step in understanding the FT-LDOS is to examine 
the symmetries of p(q, cu) by utilizing the symmetries of 
H(k) and V. We do not a priori assume any symme- 
try for the static impurity potential V. Nevertheless, 
any V is capable of being decomposed into superposi- 
tion of eight types of static impurities. Each type is 
marked by three Z2-indices, 77 = {r]x,VyiVT) which we 
define as M^zVM-J^ = rj^V, My^VMy^ = r]yV and 
UtVU^'^ = tjtV^ . Each unique combination of Z2- 
indices represents a particular type of impurity. For ex- 
ample, T] = (1,1,1) describes a static point charge im- 
purity, invariant under both mirror refiections and TRS. 
The combination rj = (1,-1,1) is an impurity invari- 
ant under mirror refiection about xz-plane and TRS yet 
changing sign under a mirror refiection about yz-p\dine. 
The physical realization for this combination of indices 
is an electric dipole which is oriented along x-axis. From 
these symmetries, we derive the exact and approximate 
symmetries of p(q, cj) and summarize them in TableU 
(for derivation, see Supplementary Materials). Here, 'ap- 
proximate' means it is valid only when [[^(cj)!^!! is small 
enough in the weak impurity limit. 

Besides the global symmetries in p(q, uj) derived above, 
the most important features of the FT-LDOS is the pres- 
ence of singularities at certain q- vectors. We start from 
writing G(k, uj) in terms of projectors: 

Tr[G(k,k + q,cj)] = 

^ Tr[Pn{k^q)Pm{k)T{uj)] 

^{uJ- Em{k) + ir]){uj - ^n(k + q) + ir?) ' ^ ^ 

where n is the band index and Pn = \un{k)){un{k)\ is 
the projector onto the n-th band. In Eq.(|4]), we have 
ignored a term Tr[Go(k, cj)](5qo, which only contributes 
to the uniform part of DOS. The numerator on the RHS 
of Eq. (HI) depends on the wave functions at k and k + q, 
while the denominator solely depends on the band disper- 
sion. The major contribution to this quantity comes from 
states that are close to the energy contours as defined 
by £^m(k) = UJ and £^n(k + q) = cj, corresponding to k- 
points where the denominator almost vanishes. Specially, 
we expect p(q, cj) to diverge at q = q^ if at some k^i, 
we have: (i) Em{Ki) = ^n(ks2 = k^i + q) = a; and (ii) 
the tangential lines of the iso-energy contours at k^i and 
ks2 are parallel. Under these conditions, ksi^s2 are re- 
ferred to as a pair of stationary points. However, through 
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TABLE I. Different types of point impurities classified by 
whether they change signs under mirror symmetries and TRS. 
From the left: the ?7-index, the possible physical realization, 
the exact symmetries of p(q,cj), the approximate symmetries 
of p(q, (ju) to the first order of V and the corresponding 
in the four-band effective model (see text). 

Eq.dH, we notice that a singularity at qg = ^(ksi — k52) 
only exists when the numerator, Tr[Pj{ks2)Pi{ksi)T{uj)] 
is nonzero. Therefore, if the two eigenstates \um{ksi)) 
and \un{'ks2)) are orthogonal to one another, we have 
Pm(ksi)Pn(ks2) = and the corresponding singularity 
in FT-LDOS at q^ must vanish. In this way, we find that 
the absence of singularities in the FT-LDOS pattern is a 
symmetry protected property of the system, because in 
general {um{ksi)\un{'ks2)) = can only be guaranteed 
by the presence of an underlying symmetry. One may 
compare the expression of p(q, cj) in Eq.([Tl|4]) with that 
of the joint density of states at q 

pj(q,cj)= (5) 

_ V - 1. 

Ntt {u - Em{k) + ir]){u - Enik + q) + ir]y 

One can see that the only difference is that in pj(q, a;) 
the factor Tr[Pmik)T{(jj)Pn{k -\- q)] is missing. Hence 
pj(q, cj) is completely determined by the geometric fea- 
tures of the iso-energy contour: at every wave vector q 
joining two stationary points, there is an intensity peak 
in pj(q, cj). A comparison between p(q, cj) and pj(q, cj) 
thus gives information of the wave functions, because any 
stationary vector that is present in the latter but absent 
in the former indicates that the two wave functions at 
the contributing stationary points are orthogonal^ ^. 

In a system with Mxz-, Myz and TRS symmetries, 
there are three types of symmetry protected orthogonal- 
ity: TRS gives us the orthogonality between two states 
of a Kramer's pair; M^^ i^yz) ensures the orthogonal- 
ity between two states with different eigenvalues of M^^ 
{Myz). The orthogonality provided by TRS leads to the 
forbidden back scattering which is a hallmark of TRS 
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topological materials. The orthogonality ensured by mir- 
ror symmetries, however, leads to new forbidden scatter- 
ing channels. In order to identify them, we need the 
mirror eigenvalues of the surface states on an iso-energy 
contour. As an example, we obtain these values here for 
Pbi-ajSuajTe. For the mirror reflection symmetry about 
xz-plane, we have [M^z^ H{kx^O)] = [Mxz^ H{kx^7T)] = 0, 
therefore each band along TX or YM may be labeled by 
its Mxz eigenvalue, either or —i. Along the TX /YM 
direction, two bands with the opposite M^z eigenvalues 
can cross each other as a gap opening in the bands is 
forbidden by the mirror symmetry (FiglTJa)). Along 
ry/XM-direction, each band can similarly be labeled 
by its Myz eigenvalue, and along the XM direction {ky- 
direction) , two bands are allowed to cross each other only 
if they have opposite eigenvalues of Myz. The two Dirac 
nodes at X (marked by 1^3,4) are protected by TRS and 
therefore appear in both directions; the other two Dirac 
nodes, D12 are protected only by M^z symmetry if they 
have opposite eigenvalues of M^z- On the other hand, 
from the absence of Dirac nodes along XM /VY and that 
the two bands bend away from each other, they can ei- 
ther have the same eigenvalues of Myz and hence repel 
each other (Fig(T]^b)), or they can have the same eigen- 
value of Myz but opposite effective mass (FiglTJc)). We 
only know that the two bands crossing at D^ /^ must have 
opposite Mxz eigenvalues. To see the Lifshitz transition, 
we find that for £^2c < ^ < ^^ic, the bands along XM 
do not cross cj, resulting in an iso-energy contour of two 
separate ellipsoids centered at D12 respectively. From 
Figdja), along Xf , the two points on the same side of 
X have opposite eigenvalues of Mxz- While for uj > Ed 
or uj < Ec2^ the bands along XM cross uo twice, resulting 
in an iso-energy contour of two concentric loops centered 
at X. From FiglTJb, c), along XM, the two points on 
the same side of X have the same eigenvalues of Myz- 

To understand how this symmetry based QPI analysis 
applies to the FT-LDOS for Pbi-^^Sua^Te, we consider the 
iso-energy contour before the Lifshiitz transition, shown 
in Figmd). Here \u{\^i)) and |ii(k2)), as well as \u{\^^)) 
and |7i(k4)), are eigenstates of Mxz each with different 
eigenvalues forcing them to be orthogonal. Now raise the 
energy, cj, such that we pass the critical point of the Lif- 
shitz transition, as seen in Figdje), where again |^i(k'i)) 
and \u{k!^)) are orthogonal as are |'^(k2)) and \u{k!^)). 
On the other hand, along f ?, all of the states are eigen- 
states of Myz^ and \u{k'^)) and \u{k!j)) are orthogonal, 
while l^i(kg)) and \u{k!^)) are also orthogonal, because 
they each have different eigenvalues of Myz. The vanish- 
ing singularities guaranteed by the mirror symmetries in 
/)(q, cj) are: qi2, q34 and q23 before the Lifshitz transi- 
tion (qa5 defined as qa-qb), and qi/2/, q3'4', q2'3', qs'T', 
q6'7' and qs's' after the Lifshitz transition. Mark that 
qi4, q23, qi'4', q2'3', qs'S', ^.Q'7' are also forbidden by 
TRS as they are scattering vectors connecting Kramer's 
pairs in the presence of TRS. 

Thus far, we have extracted the symmetries of the sin- 
gle impurity QPI pattern and located all vanishing singu- 




FIG. 1. Schematic band dispersion of the surface bands at 
vicinity of X along (a) /c^-axis and (b) /c^-axis. Red/blue 
(lighter/darker grey) means the band has eigenvalue -\-i/—i 
of Mxz in (a) and Myz in (b,c). Iso-energy contour calculated 
for the model described in the text at cj = — O.OleV in (d) and 
UJ — — O.lleV in (e). Eigenvalues of Mxz at points 1, 2, 3, 4 are 
marked beside each point, and eigenvalues of Myz at points 
5^ 6^ 7^ 8^ are similarly marked. The eigenvalues of Mxz at 
points 1^ 2^, 3^ 4^ are identical to those at 1, 2, 3, 4. 



larities guaranteed by the presence of mirror symmetries, 
as one expects in Pbi-a^Sn^^Te. However, in order to be 
more concrete, we develop a k • p- model which qualita- 
tively captures the surface band dispersion both before 
and after the Lifshitz transition and calculate the p(q, uj) 
around a point impurity. Near the Dirac points I^i,2, 
the system can be approximated by two separate two- 
band k • p models. Nevertheless, in order to describe the 
Lifshitz transition where the two cones cross each other, 
one needs to consider the hybridization between the two 
cones and, thus, the two band model must expand to four 
as the minimal number of bands in the effective Hamil- 
tonian. At X there are two degenerate states at D^, with 
El > Ed and two other degenerate states at D/^ with 
E2 < E]j. Within each doublet, one is an eigenstate of 
Mxz with an eigenvalue of +z while the other —i. We 
thus select our basis to be: (1,0,0,0) and (0,1,0,0)^ 
with energy Ei and mirror eigenvalues of Myz = +i and 
Myz = -i respectively; (0,0,1,0)^ and (0,0,0,1)^ with 
energy E2 and similarly corresponding mirror eigenval- 
ues of Myz = -\-i and Myz = —i respectively. There is 
a gauge freedom in choosing the matrices for Mxz and 
Ut under this basis, and we choose Mxz = ^^lo and 
Ut = ^^20, where = CFa ^ (^(3- Applying symmetry 
constraints, we find that, to the linear order in k, E30 
and Sio couple to the zeroth order, Sqi, 11,31 couple to 
kx and So3, 13,33 couple to ky. The effective Hamiltonian 
can hence be written as 

i^(k) = mE3o + m'Eio + (^la^^^oi + ^2x^11 + v:ix^zi)kx 

+ (^12/^03 + ^22y5]i3 + Vsy^ss)ky. (6) 
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Qualitatively fitting the data for Pbi-a^Sna^Te given in 
Ref.jil, we choose the parameters in our Hamiltonian 
to be: {m,m'} = {-0.11, -O.ljeV, {vix,V2x,vsy} = 
{ — 1.5, — 0.5, — 3}eV-A and {vsx^'^iy^'^2y} 0- Here we 
define ~ as a parameter that is very small but nonzero 
in order to avoid accidental symmetries. We represent V 
for impurities with different 7^-indices in a similar fashion 
using this basis, given in the last column of Table |T1 

Using the effective Hamiltonian for Pbi-a^Sn^^Te, we 
calculate p(q, cj) with V = O.leV * Eqo, for a static 
point charge impurity, as this is the most likely impu- 
rity to be present in Pbi-a^Sn^^Te. In Figl2fb,d), we 
plot |VV(q,^)l at a; = -O.OleV" and uj = -O.lleV re- 
spectively. (|V^p(q, cj)| is plotted in lieu of p(q, cj) since 
it shows the singularities more clearly.) In comparison, 
we make the same plots for the joint density of states, 
I V^pj(q, cj)|, in Figl2fa, c). In Figl2fa, c), we mark sev- 
eral qs's that join stationary k-points marked in FigdJ 
and at all these vectors, there are peaks in | V^pj(q, cj)|. 
However, in Figl2fb, d), we see that all peaks that are 
forbidden by either TRS or mirror symmetries vanish, 
confirming the previous analysis based on mirror eigen- 
values of the k-points. 
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protected features of the QPI. The surface is considerably 
doped with Pb atoms then, if the scattering off of surface 
electrons by these atoms cannot be ignored, the QPI pat- 
tern is generated by many similar impurities instead of a 
single impurity. However, given the random distribution 
of impurities, one can measure the autocorrelation func- 
tion of i?(r, a;), i.e., C(r,cj) = (/ dr'^ R{t^t\uj)R{t' ,uj)), 
where R{t^uj) is the FT-LDOS contributed by all im- 
purities and means the ensemble average over impu- 
rity configurations. As the Fourier transform of C{r^uj) 
is simply proportional to |p(q, to the lowest order 
of V (See Supplementary Materials), the QPI structure 
around a single impurity is refiected in the autocorrela- 
tion of the overall QPI pattern in the limit of weak and 
identical random impurities. 

It is necessary to point out that we have assumed that 
p(q, cj) contains contributions only from scattering be- 
tween states near X. However, the full QPI pattern also 
consists of the scattering among states near Y and the 
inter-pocket scattering between states near X and those 
near Y. Yet due to C4-symmetry, the contribution to 
/)(q, cj) made by states near Y can be obtained easily 
by making a fourfold rotation of p(q, uj) contributed by 
states around X. The inter-scattering between states 
near Y and those near X only contributes to p(q, uj) at 
q ~ X — Y" = (tt, tt). As these states are sufficiently sep- 
arated in q-space from features of in the small q-region, 
they are be neglected in this work. 

In samples with x ~ 1 or pure SnTe, a spontaneous 
lattice distortion break one of the mirror symmetries (as- 
sumed to be Myz) and C4 symmetry. Assuming that 
the distortion breaks M^z, we know that: first, the two 
Dirac nodes along TX are no longer protected and open 
a gap of 2E]j] second, the Lifshitz transition happens 
at two different energies, Ecx and Ecy, for the surface 
states around X and around Y. As energy decreases 
from Ef^ there are four phases: (i) no FS around X and 
two ellipsoids around X when \uj\ < Ejj^ (ii) two ellip- 
soids of FS around X and two other around Y when 
Ed < \uj\ < \Ecy\, (iii) two ellipsoids around X and two 
concentric loops around Y if \Ecy\ < |cj| < \Ecx\ and (iv) 
two concentric loops around X and two others around Y 
if |cj| > \Ecx' All these phases can be identified by QPI 
(see Supplementary Materials). In the pattern, all for- 
bidden scattering vectors forbidden by Mxz are present 
while those forbidden by Myz or TRS remain absent. 



FIG. 2. The absolute value of the divergence of the joint 
density of states, | V^/>j(q, cj)!, is plotted for cj = — O.OleV 
(a) and —O.lleV (c). The absolute value of the divergence of 
FT-LDOS around a static impurity potential V = 0.1ey*Soo, 
|VV(q,^)|, is plotted for oj = -O.OleV (b) and -O.lleV (d). 

p(q, cj) as measured in Pbi-a^Sn^^Te may contain addi- 
tional complications in identifying the crucial symmetry 
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I. DERIVATION OF THE GLOBAL SYMMETRIES OF p(q,a;) FOR DIFFERENT TYPES OF SINGLE 

IMPURITIES 

In this section we explicitly derive the symmetries of p(q, a;), exact and approximate, listed in Table I of the main 
text for single impurities with eight possible 77-indices. 

From Eq.(l), p(q, cj) is determined by the retarded Green's function G(k, k + q, a;); yet from Eq.(2), G(k, k + q, a;) 
is given in terms of the free particle Green's function Go(k, cj) and the T-matrix, T{uj). The former has exactly the 
same symmetries as the homogenous Hamiltonian: 

M^,Go{k^,ky,uj)M-^^ =Go{k^,-ky,uj), (7) 
MyzGo{kx,ky,uj)My^^ = Go{-kx,ky,uj), 
UtGo{Kuj)U-^ =G^{-k,uj). 

The symmetries of T{uo)^ however, depends on symmetries of both the system and the impurity, i.e., the 7^-index. 
Using definitions of r^-index and Eq.(3), one obtains: 

M,,T{uj)M-^ = r],V{I - Vxg{oo)V)-' = T{uj) if 7?, = 1, (8) 

^ -T{uj) if r]x -1, 

My,T{uj)M-' = rjyV{I - VygHV)-' = T{uj) if rjy = 1, 

^ -T{uj) if r]y = -1, 
UtT{uj)U-' = {r]TV{I - r]Tg{o^)V)-Y = T^{uj) if t^t = 1, 

^ —T^iuj) if TjT = —1. 

Here '~' means that the relation becomes exact only if T{uo) is replaced by its leading order, i.e., T = U. 
Combining equations in Eqs.(j8]), one easily obtains other two equations: 

{M,zMy,)T{uj){M,,My,)-' = VxVyV{I - VxVy9{^)V)-' = T{uj) if 7/,% = 1, (9) 

-T{uj) if rjxVy = -1, 

{M,,My,UT)T{uj){M,,My,UT)-^ = [r],r]yr]TV{I - r],r]yr]Tg{uj)V)-Y = '^^VxVyVr = 1 (10) 

- -T^{uj) if VxVyVT = -1. 
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Now we use Eqs.Q and Eqs.(j8]) and have 

p{q^,qy,uj) = ^ ^[Tr{M^^Go{k^,ky,Ld)T{Ld)Go{k^ ^ q^^.ky ^ qy,uj)M~^) - (q ^ -q)*] (11) 

k 

= ^[Tr(Go(fc^, fcy,u;)M,,T(u;)M-/Go(fc, + q,, ky - qy,w)) - (q ^ -q)*] 

k 

= p{Qx,-qy,^) if = 1 

^ -p{Qx, -Qy,^) if Vx -1, 
p{q:,,qy,u) = ^[Tr{My,Go{k^,ky,uj)T{uj)Go{k^ ^ q^.ky ^ qy,uj)M-^) - (q ^ -q)*] (12) 

k 

= -^^[Tr{Go{k^,ky,u)My,T{u)M-^^ - (q ^ -q)*] 

k 

= p{-qx^qy^^) if% = 1 
~ -p{-Qx^qy^^) if Vy = -1, 

p{qa,,qy,uj) = V'[Tr(/7TGo(/Ca^, A;^,cj)T(cj)Go(A;^ q^^.ky ^ qy,uj)U~^) - (q ^ -q)*] (13) 

k 

= ^[Tr(Go(-fc^ - q^,-ky - qy,w){UTT{Lo)U-^Y Go{-k^, -ky,uj)) - (q ^ -q)*] 

k 

= -^J2i^r{Go{k, - q,,ky - qy,Lo){UTT{co)U-TGo{k,,ky,oj)) - (q ^ -q)*] 

k 

= piQx^Qy^^) if r/T = 1 
~ -p{Qx,qy,^) if = -1. 

We can also use Eqs.© and Eqs.(|9]) to obtain, very similarly, 

p{Qx,qy,^) = p{-qx,-qy,^^) VxVy 1 (14) 

^ -p{-Qx,-qy,^) ^^VxVy = -1, 
p{qx,qy,uj) = p{-qx, -Qy,^) if VxVyVr = 1 (15) 

^ -p{-Qx, -Qy,^) if ^x%^T = -1. 

Up to now we have derived all symmetries, exact and approximate, of p(q, uj) listed in Table I. 

II. AUTOCORRELATION OF R{r,uj) IN THE WEAK IMPURITY LIMIT 

In the text, we mention that although we have only studied the case of a single impurity, the result extends to 
the case where there are random impurities of the same type and strength and when the impurity strength is weak. 
Here we prove that under these conditions, the FT-LDOS around a single impurity, p(q, cj), can be related to the 
autocorrelation function of the overall LDOS resulted from many impurities, R{t,uj). By definition, 

R{r,cu) = V \M^)\^^{co - En) = --/m[G(r,r,c^)], (16) 

Z ^ 

n 

where V^n(r) is an eigenfunction of the Hamiltonian, having eigenvalue En- The Fourier transform of R{t) is, in terms 
of the retarded Green's function, 

i?(q, uj) = ^Yl Tr[G{K k + q, a;) - G* (k + q, k, uj)] , (17) 



k 



which is identical to the expression of p{r) in the case of a single impurity. But now G(k, k + q, cj) does not have a 
simple expression as in Eq.(3) in the text, since the system contains more than one impurity. 
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By definition, we tiave the Fourier transform of the correlation function C(r) as (from here, we suppress the argument 
of uo for concision) 

C(q) = I d'r'dMRir + r')Rir') exp(-iq • r)) (18) 

= (l^?(q)l')- 

Therefore the Fourier transform of the autocorrelation function C(r) is exactly the ensemble average of the module 
squared Fourier component R{q). 

For an arbitrary configuration of identical impurities, the retarded Green's function is 

G(k, k + q) = Go(k)(5qo + Go(k)F(q)Go(k + q) + ^ ^ Go{k)V{k - k')Go{k')V{k' - k - q)Go(k + q) + ...(19) 

where T^(q) = F^-exp(— iq • r^) is the Fourier transform of the impurity potential. From the above equation we 
have i^(q) as 

i?(q) = Po5^o + 4? E Tr[Go(k)F(q)Go(k + q) - G5(k + q)F(q)G5(k)] (20) 



k 



^2 E E Tr[Go{k)V{^ - k')Go{k')V{k' - k - q)Go(k + q) - GS(k + q)F(k' - k - q)GS(kOF(k - k')GUk)] 



k,k' k 



For q 7^ 0, the lowest order in |p(q)P is the second order in V: 

|i?(q)P = |Tr{l E[Go(k)^(q)Go(k + q) - G5(k + q)F(q)GS(k)]}p| Eexp(-iq • r.)|^ (21) 

k i 

= Kq)n^exp(-iq.r,)|^ 



Finally, we use the independence (randomness) of the distribution of to prove that 

(I ^ exp(-zq . rOl') = A^i + ^(exp(zq • (r, - r,))) = AT,, 

where Ni is the total number of impurities. Up to this point, we have proved that Fourier transform of the auto- 
correlation function C{r) of the overall LDOS R{t) is proportional to |p(q)P in the limit of weak impurities. In 
other words, the singularities in p(q) can be revealed in the autocorrelation of LDOS in the case of multiple identical 
impurities. 



III. QPI WITH RHOMBOHEDRAL DISTORTION 



In the text, we have mentioned that QPI can be used to identify the four types of iso-energy contours at different 
energy u. In this section we use the effective model as example to see how p(q, cj) evolves as u varies. 

Without loss of generality, we assume that the distortion breaks M^^ but preserves Myz TRS. In the presence of 
the rhombohedral distortion, the fourfold rotation symmetry is broken and we have to consider the contribution from 
states near X and Y separately. Any term that breaks M^^ but preserves Myz and TRS is of index 77 = (—1, +!,+!). 
According to Table I, in the effective theory around X, this term is AS23. Without distortion, the effective theory 
around Y is obtained by rotating the model in Eq.(6) by 7r/2, and in this theory, the distortion term corresponds 
to AH23. In Figini we choose A = O.leV and calculate the iso-energy contours at = — O.OleV in (a), — O.OTeV in 
(b), — O.lleV in (c) and — 0.15eV in (d). From these figures, the change in the topology of the contours is clear: it 
changes from only two separate ellipsoids around Y in (a), to two separate ellipsoids around both X and Y in (b), to 
two separate ellipsoids around X and two concentric loops in (c), and finally to two concentric loops around both X 
and Y at highest u in (d). The resultant QPI patterns, |V^p(q, a;)|, at these energies are plotted in FiglSl^e, f, g, h) 
respectively. These patterns are simply superpositions of patterns contributed by states around X and those around 
F; each pattern explicitly breaks fourfold rotational symmetry. Another important feature is that all q- vectors that 
are forbidden by Mxz are now present. For example, in FiglSl^f), the inner circle is contributed by the surface states 
around X. On this circle, the left and right ends, defined as q^ and q^^ are due to the scattering between k's at 
the two ends of either one of the two pockets near X. When M^z is unbroken, these two k's have opposite Mxz 
eigenvalues and hence scattering between them is forbidden. But with Mxz broken, qL,i? are allowed by symmetry. 
Compare this with Fig. 1(b), we can see that the distortion makes the intensity at q_L,R finite. 
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FIG. 3. (a,b,c,d): The iso-energy contours at cj = —O.OleV, — O.OTeV, —O.lleV and — 0.15eV, respectively, in the phase with 
rhombohedral distortion, breaking Mxz symmetry. (e,f,g,h): V^/>(q, cj) calculated for he surface states in the rhombohedral 
phase around a single charge impurity at the same energies listed above. 



